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Abstract Slow-light enhanced optical detection in liquid-infiltrated pho- 
tonic crystals is theoretically studied. Using a scattering-matrix approach 
and the Wigner-Smith delay time concept, we show that optical absorbance 
benefits both from slow-light phenomena as well as a high filling factor of 
the energy residing in the liquid. Utilizing strongly dispersive photonic crys- 
tal structures, we numerically demonstrate how liquid-infiltrated photonic 
crystals facilitate enhanced light-matter interactions, by potentially up to 
an order of magnitude. The proposed concept provides strong opportunities 
for improving existing miniaturized absorbance cells for optical detection in 
lab-on-a-chip systems. 



1 Introduction 



Optical techniques are finding widespread use in chemical and bio-chemical 
analysis, such as absorbance, fluoresce nce, Raman-scattering, a nd surface- 
plasiri o n-resonance n i easure r nents, see Moeenseri et al. (2QQJ); Hess et al 
(|2002l ):'Kneip p et all (|2002h : iHomola et all (|l999f l. In particular, the Beer- 



Lambert absorbance measurement has become one o f the cl assical quantita- 
tive workhorses in analytical chemistry, ISkoog et all (|l997l) . The increasing 
emphas is on miniatur ization of chemical analysis systems during the past 
decade, 



Janasek et al 



(|2006) . has naturally stimulated a considerab l e effor t 
in integrating microfluidics. Squires and Ouake J 2005 )^ Whitesides ( 20061). 
and optics in lab-on-a-chip microsystems, Verpoortd (|2003 ): Mortensen et al. 
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Fig. 1 Panel (a) illustrates a typical macroscopic cuvette while panel (b) shows 
a microscope image (top-view) of an equivalent lab-on-a-chip implementation of 
a microfluidic channel (vertical direction) integrated with planar optical waveg- 
uides (horizontal direction). Courtesy of K. B. Mogensen and J. P. Kutter (MIC 
- Department of Micro and Nanotechnology, Technical University of Denmark, 
www.mic.dtu.dk) . 



(j2007l ). and more recen tly this has to a lar g e extent powered th e emerging 



field of optofluidics, see lPsaltis et al. I (|2006l ): lMonat et all (|2007[) . However, 



lab-on-a-chip implementations of the above mentioned optical techniques 
call for new approaches going beyond a simple miniaturization of existing 
optical setups. The need for a new perspective is perhaps best appreciated 
by emphasizing the optical drawback of miniaturization. Reducing the opti- 
cal path length will eventually also decrease the sensitivity significantly. In 
other words, the sensitivity fundamentally scales with volume hosting the 
interaction between light and the solution of analytes. 

Panel (b) in Fig. [1] illustrates a typical lab-on-a-chip implementation of 
an absorbance cell. The optical path length L is often reduced by several 
orders of magnitude compared to typical macroscopic counterparts, such as 
the cuvette shown in panel (a). A typical size reduction by two orders of 
magnitu de will penalize the op tical sensitivity in an inversely proportional 



manner, iMogensen et al. I (|2003h 



In the above context, photonic crystals (PhC) are potential candidates 
for the desi red new approach on th e light-matter interactions as recently 
reviewed bv lMortensen et al ( 2007 ). The concept of PhCs refers to a class 
of artificial electromagnetic str uctures offering high l y engineered dispersion 
properties, as first suggested by I Yablonovitch ( 1987t ): [johnl ( 19871 ). PhCs are 
porous structures, but with the voids arranged in a highly regular way. Thus, 
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for a wavelength A comparable to the periodicity A of the PhC, the elec- 
tromagnetic radiation interacts strongly with the matter. PhCs have over 
the years been recognized as strongly dispersive environments supporting a 
number of p henomena including photo n ic band g aps a nd slow-light prop- 
agation, see Ijoannopo ulos et al. ( 19951 ); Sakoda ( 2005 ). The key-point in 



the present context is that liquid-infiltrated PhCs will overall inherit the 
unusual dispersion properties from their non- infiltrated counterparts, thus 
also changing the way light interacts with bio-molecules dissolved in the liq- 



iig- 

uid. iMortensen et al.l ([2007). The effects can be quite pronounced compared 



to light-matter interactions in a spatially homogeneous liquid environment 
and in this paper we follow this route to compensate the reduced optical 
path. 

In particular, we consider slow-light enhanced optical detection in liquid- 
infiltrated photonic crystals, thus potentially compensating for the cost of 
miniaturization and reduction in optical path length. 

The enhancement factor 7 = a/ai for an absorbing liquid (with absorp- 
tion coefficient ai) can be expressed as 

r c/ni 

7 = /x , (1) 



as recently derived by ^Mortensen and Xiaol ( 2007 ). Here, < / < 1 is 



a dimensionless number quantifying the relative optical overlap with the 
liquid. The fraction on the right-hand side expresses the ratio of the group 
velocity c/n; in the bare liquid to the group velocity Vg in liquid-infiltrated 
photonic crystal, thus clearly illustrating the enhancement by slow-light 
propagation (vg <C c). In the following, we first derive the above expression 
before we with numerical simulations illustrate the slow-light enhancement 
in liquid-infiltrated PhCs with three different periodic structures. 



2 Theory 

2.1 Electromagnetic wave equation 

The electromagnetic problem that we consider is governed by the electro- 
magnetic wave equation for the electrical field, 

2 

V X V X 1^;) = £^^|£;), (2) 

where for the absorbing liquid the dielectric function e has a small imaginary 
part, uj is the a ngular frequency, and c is th e velocity of light in vacuum. In 
a recent paper, IMortensen and Xiaol (|2007h derived Eq. ^ rigorously with 



the aid of standard first-order electromagnetic perturbation theory 

Lu(E\6e\E) 

5^ = -77 \ I I / (3) 
2 {E\e\E) ^ ' 
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and m ore details may be found in our recent work, see iMortensen et alJ 
(j2007n . In Eq. ([3]), Suj is the first-order frequency shift introduced by a 
perturbation 6e of the dielectric function e (assumed frequency independent) 
with the electric field being the unperturbed one. In the present context the 
perturbation is associated with absorption in the liquid-part of space thus 
giving rise to an imaginary frequency shift. In the following we emphasize 
this by explicitly writing the shift as iSuj with Suj being real. For the temporal 
dependence of the intensity we have I{t) = /o|e~"^*| = p-^Sijt we 
may define an absorption rate given hy F = 2Suj. 



2.2 Scattering-matrix approach 



Here, we will emphasize the light-matter interaction time and offer an alter- 
native derivation based on scattering matr ices combined with the concept 
of the Wigner-Smith delay time. Following iBeenakker and Brouwen (j200ll ) 
we start from the scattering matrix 



S{uj) 



Sii{uj) 5*12(0;) 

52l(w) S'22(t^) 



r{uj) t'{uj) 
t{uj) r'{uj) 



(4) 



where Sij are the S-parameters connecting incident and out-going plane- 
wave (or Bloch) states to the left and the right of the scattering region. 
Next, let So{lo) denote the scattering matrix in the absence of absorption so 
that for a weak absorption we have a small imaginary shift iSuj in frequency 
and thus S{uj) ~ 5'o(w -I- iSuj). Taylor expanding the right-hand side we get 



Sito) ~ So{lo) + iSio^^^ = Soiuj) 

OUJ 



1 - ^QoH 



(5) 



In the second equality F — 2Suj is the absorption rate and we have further- 
more used the unitarity of the scattering matrix, SqSq = 1, to introduce 
the Hermitian Wigner-Smith delay time matrix 



?o(w) = -tS^^iuj)—^—. 



(6) 



Now, we would like to calculate the transmission probability T through a 

single eigenmode in which case S* is a two-by-two matrix. Assuming time- 

I 1 2 I 1 2 

reversal symmetry we have T = \Si2\ = |'5'2i| which we have the freedom 
to rewrite as 

T{lo)^^Tt{S{lo)SHlu)}+0{1-T). 



(7) 



For near-resonance transmission, i.e. for a low reflection, we get 

2 



1 1 ^ 

T{u;) -Tr [1 - FQoicj)] = i,J2^^~ ^*"(^)] = ^ - (8) 
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where we have used the unitarity of Sq together with the cychc invariance 
of the trace. In the second equahty we have introduced the eigenvalues 
ti = t2 = T of the Wigner-Smith delay time matrix where t is often referred 
to as the dwell time. We may use perturbation theory, Eq. ([3]), to calculate 
the absorption rate which gives 

(E\D), 

r = fxri, f= \ V (9) 
{e\d) ^ ' 

where we have introduced the displacement field |_D) = e\E) and /] is the 
absorption rate of the liquid itself. For the filling factor we have < / < 1 
and the integral in the numerator is restricted to the region containing 
the absorbing fluid while the integral in the denominator is spatially un- 
restricted. The presence of the filling factor / is quite intuitive since only 
the fraction / of the light residing in the fluid can be subject to absorp- 
tion. Note that co mpared to the definition no rmally used by the photonic 
crystal community, Joannopoulos et al. ( 1995[) . our definition quantifies the 



relative overlap with the liquid-part of space rather than the dielectric 
part given by 1 — /. Finally, we rewrite Eq. ([8|) by re-exponentiating, i.e. 
(I — X) — > exp(— x), so that 

T = cxp(-rT) EE exp(-7 X FiTi) (10) 

corresponding to the expected exponential damping of the resonance due 
to absorption. In the second equality we have introduced 7 = ajai = 
[Ft) / {FiTi) and using that t = L/vg for the periodic medium and r; = 
L/{c/ni) in the bare liquid we finally arrive at Eq. p]). 

The expression for 7 clearly demonstrates how optical absorbance bene- 
fits from slow-light phenomena. The effective enhancement of the absorbance 
may thus naturally be interpreted as an effective enhancement of the light- 
matter interaction time. These conclusions may be extended to also non- 
periodic systems, including enhanced absorbance in disordered sys t ems a s 
well as intra-cavity absorbance configurations, see iMortensen et al. 1 (|2003). 



In general, fluids/liquids have a smaller refractive index than the typical 
solid materials comprising the photonic crystal itself. Thus, increasing the 
filling factor (and thereby the light-liquid interaction) implies the confine- 
ment or guidance of light in the low index regions, which is typically not the 
case explored in more ordinary applications of PhCs for e.g. optical commu- 
nication purposes. However, a high light-liquid overlap can be achieved in 
some structures, e.g., slot waveguides, hollow anti-resonant reflective opti- 
cal waveguides, as well as photonic-crystal structures. For liquid-infiltrated 
photonic crystals and photonic crystal waveguides, it is possible to achieve 
Vg <^ c and at the same time have a filling factor of the order unity, / ~ 1, 
whereby significant enhancement factors become feasible. 
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Fig. 2 (a) Photonic band structure for normal incidence of either TE or TM 
polarized light on a Bragg stack of period yl = aj + 02 with n; = 1.33, 712 = 3, 
ai — 0.6863/1, and 02 = 0.3137/1. Photonic band gaps are indicated by yellow shad- 
ing and the dashed line indicates the long- wavelength asymptotic limit, Eq. (|lip . 
(b) Corresponding enhancement factor which peaks and exceeds unity close to the 
photonic band-gap edges. The dashed line indicates the long-wavelength asymp- 
totic limit, Eq. (fT^ . 



3 Numerical examples 

Having derived a formal expression for th e enhancement factor 7, we next 
apply a full- vectorial plane wave method of Johnson and Joannopouloi ( 2001 ) 
to calculate the eigenmodes of Eq. ^ in the non-absorbing limit and from 
these solutions we evaluate the enhancement factor, Eq. ([1]), for a number 
of geometries. 



3.1 Bragg stack 

Let us first illustrate the slow-light enhancement for the simplest possible 
structure which is a Bragg stack with normal incidence of electromagnetic 
radiation, as shown in the inset of panel (b) of Fig. 2. The structure is 
composed by the low-index material layers of width a; — 0.6863/1 being a 
liquid with refractive index ni = 1.33 and the high-index layers have a width 
02 = 0.3137/1 and a refractive index 712 = 3. Panel (a) of Fig. 2 shows the 
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Fig. 3 (a) Photonic band structure for propagation of TM polarized light along 
the T^X direction in a square lattice of period A with dielectric rods of diameter 
d = 0.2906/1, n = 3.24 and n; = 1.33. Photonic band gaps are indicated by yellow 
shading, (b) Corresponding enhancement factor which exceeds unity for the flat 
bands in general and the fourth band in particular. The right y-axis shows the 
results in terms of the free-space wavelength when results are scaled to a structure 
with A = 1000 nm. 



corresponding band structures, where photonic band gaps are indicated by 
yellow shading and the dashed line indicates the long-wavelength asymptotic 
limit to a close-to-linear dispersion 



ckA 



aim + 02^2 



(11) 



When approaching the band-gap edges, the disp ersion flattens. Note tha t 
the first band is a dielec tric-like one, see e.g. iMortensen et al.l (|2007l ). 



Joannopoulos et al.l (|l995r ). meaning that most of the energy localizes in 
the high dielectric region. For the first band it is well-known that the flat 
dispersion originates from a spatial localization of the field onto the high- 
index layers and thus / ^ 1 near the band edges where the inverse group 
velocity diverges. However, in spite of the localization, the enhancement 
factor may still exceed unity as shown in panel (b) where the dashed line 
indicates the long-wavelength asymptotic limit with 



/ 



(1 
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Fig. 4 (a) Photonic band structure for propagation of TE polarized light along 
the rX direction in a PhC waveguide in a triangular lattice of period A with holes 
of radius r — 0.37/1 and n — 3.24. The PhC waveguide is formed by enlarging 
the radius (ri = 0.50/1) of the air holes in a single row, which are infiltrated by 
a liquid with n; — 1.33. The grey shading indicates the regions of finite density 
of states in the surrounding PhC while the yellow shading shows the waveguide 
mode-gap region, (b) Corresponding enhancement factor 7. The inset shows the 
structure of the photonic crystal waveguide. 

For the second band, the enhancement factor still exceeds unity, which is 
attributed both to slow-light, as well as to the large filling factor. The filling 
factor is about 20% at the band edge for the second band. In order to further 
benefit from the slow-light enhanced light-matter interaction, we obviously 
have to pursue optofiuidic structures supporting both low group velocity 
and at the same time large filling factors. 

3.2 Pillar-type photonic crystal 

Figure 3 shows the case for a 2D square photonic crystal with high-index 
dielectric rods in liquid. Compared to the Bragg stack in Fig. [21 some of 
the modes in this structure have both a low group velocity and at the same 
time a reasonable value of the filling factor /. In particular, the fourth band 
in panel (a) is quite flat. The corresponding enhancement factor 7 exceeds 
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4 for a bandwidth of order 100 nm for a pitch around A ^ 1000 nm, as 
indicated on the right y-axis in panel (b). 



3.3 Void-type photonic crystal waveguide 



Recently, it was shown bv lGersen et al. (|2005h how photonic crystal waveg- 



uides can be used to slow down the propagation velocity of light. Here, we 
illustrate how this may be utilized for slow-light enhancement in a PhC 
waveguide structure. We consider the PhC waveguide shown in the inset 
in panel (b) of Fig. IH where the PhC is two-dimensional with a triangular 
array of holes in an InP/GalnAsP dielectric background with a refractive 
index of n = 3.24. The radius of the air holes is r = 0.37yl. In order to 
realize a non-solid waveguide, the PhC waveguide is formed by slightly en- 
larging the radius (ri ) of the air holes in a single row which are subsequently 
utilized for liquid-infiltration. In order to obtain a high filling factor, we op- 
timize the waveguide structure and obtain / « 40% when ri = 0.5/1. The 
enhancement factor exceeds unity near the band edge, which is attributed 
to both slow-light propagation as well as to the large filling factor. 



4 Conclusion 

In conclusion, we have studied the potential of using liquid-infiltrated pho- 
tonic crystals to enhance optical absorbance. Using a scattering-matrix ap- 
proach we have re-derived a general expression for the enhancement factor, 
thus confirming the interaction-time interpretation offered in previous work. 
Our results clearly demonstrate how optical absorbance may benefit from 
slow-light phenomena combined with a filling factor close to unity. With 
the aid of photonic crystals, it is possible to design an optical enhanced 
detection system which has both a slow group velocity and a high filling 
factor. The slow-light enhancement of the absorption, by possibly an order 
of magnitude, may be traded for yet smaller miniaturized systems or for 
increased sensitivity of existing devices. 
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